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Abstract 

The branching ratios and CP violations of the B — > Dtt decays, including both the color-allowed 
and the color-suppressed modes, are investigated in detail within QCD framework by considering all 
diagrams which lead to three effective currents of two quarks. An intrinsic mass scale as a dynamical 
gluon mass is introduced to treat the infrared divergence caused by the soft collinear approximation in 
the endpoint regions, and the Cutkosky rule is adopted to deal with a physical-region singularity of the 
on mass-shell quark propagators. When the dynamical gluon mass fig is regarded as a universal scale, 
it is extracted to be around fig = 440 MeV from one of the well-measured B Dtt decay modes. 
The resulting predictions for all branching ratios are in agreement with the current experimental 
measurements. As these decays have no penguin contributions, there are no direct CP asymmetries. 
Due to interference between the Cabibbo-suppressed and the Cabibbo-favored amplitudes, mixing- 
induced CP violations are predicted in the B D^tt^ decays to be consistent with the experimental 
data at 1-cr level. More precise measurements will be helpful to extract weak angle 2(3 + -y. 

PACS numbers: 13.25.Hw, ll.30.Er. 12.38.Bx. 
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I. INTRODUCTION 



Nonleptonic 5-meson decays are of crucial importance to deepen our insights into the flavor 
structure of the Standard Model (SM), the origin of CP violation, and the dynamics of hadronic 
decays, as well as to search for any signals of new physics beyond the SM. However, due to 
the non-perturbative strong interactions involved in these decays, the task is hampered by 
the computation of matrix elements between the initial and the final hadron states. In order 
to deal with these complicated matrix elements reliably, several novel methods based on the 
naive factorization approach (FA) such as the QCD factorization approach (QCDF) 2], the 
perturbation QCD method (pQCD) [3|, and the soft-coUinear effective theory (SCET) J], have 
been developed in the past few years. These methods have been used widely to analyze the 
hadronic i?-meson decays, while they have very different understandings for the mechanism of 
those decays, especially for the case of heavy-light final states, such as the B Dn decays. 
Presently, all these methods can give good predictions for the color allowed b'^ D^tx^ mode, 
but for the color suppressed B^ D^tt^ mode, the QCDF and the SCET methods could not 
work well, and the pQCD approach seems leading to a reasonable result in comparison with 
the experimental data. In this situation, it is interesting to study various approaches and find 
out a reliable approach. 

As the mesons are regarded as quark and anti-quark bound states, the nonleptonic two 
body meson decays concern three quark-antiquark pairs. It is then natural to investigate the 
nonleptonic two body meson decays within the QCD framework by considering all Feynman 
diagrams which lead to three effective currents of two quarks. In our considerations, beyond 
these sophisticated pQCD, QCDF and SCET, we shall try to find out another simple reliable 
QCD approach to understand the nonleptonic two body decays. In this note, we are focusing 
on evaluating the B Dn decays. 

The paper is organized as follows. In Sect. II, we first analyze the relevant Feynman 
diagrams and then outline the necessary ingredients for evaluating the branching ratios and 
CP asymmetries oi B ^ Dn decays. In Sect. Ill, we list amplitudes of -B ^ Dn decays. The 
approaches for dealing with the physical-region singularities of gluon and quark propagators are 
given in Sect. IV. Finally, we discuss the branching ratios and the CP asymmetries for those 
decay modes and give conclusions in Sects. V and VI, respectively. The detail calculations of 
amplitudes for these decay modes are given in the Appendix. 
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II. BASIC CONSIDERATIONS FOR EVALUATING B ^ Dn DECAYS 



We start from the four-quark effective operators in the effective weak Hamiltonian, and then 
calculate all the Feynman diagrams which lead to effective six-quark interactions. The effective 
Hamiltonian for B Dn decays can be expressed as 

Heff = ^v,tv:,[CMO,{f,) + C2{^I)02{^^)] + h.c, (i) 

where Ci and C2 are the Wilson coefficients which have been evaluated at next-to-leading 
order j^, Oi and O2 are the tree operators arising from the ly-boson exchanges with 

Oi = icibi)v-AidjUj)v-A, O2 = {cibj)v--A{djUi)v-A, (2) 

where i and j are the SU(3) color indices. 

Based on the effective Hamiltonian in Eq. ([1]), we can then calculate the decay amplitudes 
for B D+Ti", B L)°7r°, and B~ D^tx~ decays, which are the color- allowed, the 
color-suppressed, and the color-allowed plus color-suppressed modes, respectively. All the six- 
quark Feynman diagrams that contribute to 5^ ^ D^ir^ and B^ D^n^ decays are shown 
in Figs. [T]l3] via one gluon exchange. As for the process B~ D^tt" , it doesn't involve 
the annihilation diagrams and the related Feynman diagrams are the sum of Figs. [T] and [2j 
Based on the Isospin symmetry argument, the decay amplitude of this mode can be written 
as A{B- L)°7r-) = A(B° D^n') - V2A(B° D^n^). The explicit expressions for the 
amplitudes of these decay modes are given in detail in next section. 

The decay amplitudes of i? ^ Dtt decay modes are quite different. For the color- allowed 
B^ D^7r~ mode, it is expected that the decay amplitude is dominated by the factorizable 
contribution A/ac (from the diagrams (a) and (b) in Fig. [1]), while the nonfactorizable contri- 
bution Anonfac (from the diagrams (c) and (d) in Fig. [1]) has only a marginal impact. This 
is due to the fact that the former is proportional to the large coefficient ai = Ci + ~ 1, 
while the latter is proportional to the quite small coefficient 02 = C2 + ~ 0. In addition, 
there is an addition color-suppressed factor in the nonfactorizable contribution Anonfac- 
In contrast with the B^ D'^tx^ mode, the nonfactorizable contribution Anonfac (from (c) 
and (d) diagrams in Fig. [2]) in the 5^ D^n^ mode is proportional to the large coefficient 
Oi = Ci -|- ~ 1, and even if with an additional color-suppressed factor its contribution 
is still larger than the factorizable one Afac (from (a) and (b) diagrams in Fig. [2]) which is 
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FIG. 1: The factorizable ((a) and (b)) and nonfactorizable ((c) and (d)) diagrams contributing to the 



color-allowed — > D~^tt decay. 




(b) 



(c) 



(d) 



FIG. 2: The factorizable ((a) and (b)) and nonfactorizable ((c) and (d)) diagrams contributing to the 
color-suppressed D^tt^ decay. 



proportional to the quite small coefficient 02 = C2 + ~ 0. Thus, it is predicted that the 
decay amplitude of this mode is dominated by the nonfactorizable contribution Anonfac- As for 
the B~ D^TT^ mode, since its amplitude can be written as the sum of the ones of the above 
two modes, it is not easy to see which one should dominate the total amplitude. 

The branching ratio for B — > Dir decays can be expressed as follows in terms of the total 
decay amplitudes 

^ Dn) = \A{B ^ Dn)f (3) 

where tb is the lifetime of the B meson, and pc is the magnitude of the momentum of the 
final-state particles D and tt in the 5-meson rest frame and given by 

Pc = ~ ("^^ + ^^y] I'^B - i'f^D - m^)'^] . (4) 

As is well-known, the direct CP violation in meson decays is non-zero only if there are 
two contributing amplitudes with non-zero relative weak and strong phases. The weak-phase 
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difference usually arises from the interference between two different topological diagrams. For 
three B Dtt decays, it is seen from the Feynman diagrams in Figs. [T]l3] that there are no 
weak-phase differences, and hence no direct CP violation in all these three modes, we shall 
then consider the mixing-induced CP violation. 

As the final states D^tt~ can be produced both in the decays of 5^ meson via the Cabibbo- 
favored {h c) and in the decays of meson via the doubly Cabibbo-suppressed {b — ^ u) 
tree amplitudes. The relative weak-phase difference between these two amplitudes is —7 and, 
when combining with the B^ — B^ mixing phase, the total weak-phase difference is —(2/5 + 7) to 
all orders in the small CKM parameter A. Thus, the B^{B^) D^ir' decays can in principle 
be used to measure the weak phase 7, since the weak phase (3 has been measured with high 
precision. The time-dependent CP asymmetry of such decay modes is defined as: 

^ T(b' D+vr-(At)) - r(E° D+vr-(At)) 
T(B° ^ D+7i-{At)) +T{BO ^ D+7i-{At)) 
-a, + (ce + a^^''") COS {AniB ■ t) + afj^J' sin (Am^ ■ t) 

— Sn+T^- sin(Amt) — Cd+t,- cos(Amt), (5) 

where Am is the mass difference of the two eigenstates of Bd mesons, and Sd+t^- and Cd+t^- 
are given as 

_ 2Jm(Az)+^-) _ _ 1 - IXp+n-l^ _ 



1+ A 
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with 

Where the rephase- invariant quantities a^, a^/ and cte+e' |6| characterize the indirect, direct 
and mixing-induced CP violations respectively. As <^ 1 for neutral B system, we have 
Cd^-k- — —a^^'^ which characterizes direct CP violation. Defining g and h as the amplitudes 
of — > D'^Tx^ and B'^ — > D~^7r~ decay modes, respectively, we can further express these two 
CP asymmetries as 

Cd+tt- = |,|2 , I 12 = 2 I 1 

_ -2|^||/i|sin(2/3 + 7 + 5) _ -2^sin(2/3 + 7 + (5) 

where 2; = \h/g\, and 5 = arg|/i/(7| represents the relative strong-phase difference between the 
two amplitudes g and h. 

Similarly, we can define another two CP- violating parameters C£)-7r+ and S'£)-^+ for the 
B^{^) ^ 0^71+ decays 

_l--z' _ -2zsin(2/3 + 7-^) 

1 + -2 1 + 2;^ 

with the parameter z defined as z = \h/g\ = z. Here the amplitudes h and g are the charge 
conjugations of the amplitudes h and g. Since the magnitude of the Cabibbo-suppressed decay 
amplitude \h\ is much smaller than that of the Cabibbo-favored decay amplitude \g\, the ratio 
z should be quite small and is found to be about 0.02 in our framework. Thus, to a very good 
approximation, Cd+t^- = —Co-n+ — — 1? and the coefficients of the sine terms are given by 



5^+,- = -2z sin(2/5 + 7 + 5) = , Sd-^+ = -2z sm{2/3 + 7 - 5) = • (11) 

To compare with the current experimental data, one usually define the following two quantities, 
which are given by the combination of two CP- violating parameters Sd+w- and Sn-n+j 

a = {Sd+tt- + SD-n+)/'2 , c = {So+TT- - SD-n+)/'2, (12) 
which can provide constraints on the weak phase 2/5-1-7 and the strong phase 6. 

III. B ^ Dn DECAY AMPLITUDES 

Using the methods given in the Appendix, we can get the B Dn decay amplitudes, 
which are composed of three parts: the factorizable contribution Afac, the nonfactorizable 
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contribution Anonfac, and the annihilation contribution Aanni- The amphtude of -B ^ D^tt 
mode is found to be 



with 



A{B D+TX-) = VcbV*Mfac + Anonfac + A 



G C C 

= + Tr)Tr dx dy(j)B{x)(j)D{y) 

\I I ^^c ^^c Jo Jo 



(13) 



X 



ym% + 2mbmB — 2ym£,mB + nihrni)] 



-mjnB — 2xmB)mB — 2mcm£)) 



m. 



(14) 



where Df, = m^y — ml, Dc = —m?BXX — and fc^ = m?Bx{x — y). (j)'s are the wave functions 
of mesons. For the 5-meson wave function, we shall take the form given in [3] 

2" 



0i?(p) = iVBp'(l -p)'exp 



_1 / ptub X 
2\ub J 



(15) 



with ujb = 0.25 GeV, and Nb being a normalization constant. The D meson distribution 
amplitude is given by 

0z5(y) = 6y(l-y)[l + Cz,(l-2y)], (16) 

with the shape parameter Cd = 0.8. For the tt meson light cone wave functions, we use the 
asymptotic form as given in Refs. s], [o], Q : 



(u) = Quu, 



(u) = 1. 



(17) 



with u = 1 — u. 



G C C 

Anonfac = f B fo fnT^^sifJ') {C2 + ^) ^ / dx dy dz (j)B{x)(t)D{y)(t){z) 

y2 i\c Jo Jo Jo 

{{-X + z)m% - (x - y)mDmB) -jj^ 

+ [{-2x + y + z)m\- {x -y)mDmB)jj^U (18) 
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where Dd = m^(a: — y){x ~ z) and Du = ra\{x ~ y){x — z). 

Aanni = Jb fo fnT^asifi) -r^ / dx dy dz (poiv) 
\ 2. ^^c JO JO Jo 



X 



(C2 



Ci 



Nc 

finirric + 2zmD) 



zm\ — 2mcmr))(l){z) + fiTr{mc — 2zmD + 4mi5)0^(z) 
0^(2) \ 
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6 ; /^caA;2 



-ym^(f){z) + 2jj,^mr,{y + 1)0,1(2; 



B 



-'-^ua'^a 



((x - y)mB + mf,)m%(t){z) + /i^m^, ((2x -y - z)mB 



- 4:mb)(t>jr{z) + fi^momBiy - z)- 



a{Z) \ rriB 



6 / D.±^„ 



( 0' {A 
+ ( (x - ^)m|0(2) - Unmody - 2) h (2x - y - 2)0^(2;)) 



mi 



(19) 



m%z-m:, Dua = m%y, D^a = m%{x-y){x~z)-ml, D^a = m%{x 



where 0^(z) = -Dca - "^s- 

|/)(x — z) and /c^ = m'^^yz. The annihilation contribution is found to be much smaller than the 
ones from the factorizable and the nonfactorizable diagrams. Numerically, it is negligible. 



For the color-suppressed B D n decay, its amplitude can be written as 



(20) 



with 



A 



fac 



G C C 

^/i?/DAvra.(/i)(C2 + ^)— ^ dx dz(j)B{x) 
V2 iVc iVc Jo Jo 

-zm% + 2mbm.B)(l){z) + yU^(2zmB - mh)(j)T,{z) 



+ fi^{2{z + l)mB + mi,) 



4>'^{z)\ m 
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+ 2fi^x(f)^{z) 



m 



B 



(21) 



here Dt = rn^z — ml^ Dd = m^{x — y){x — z) and /c^ = m'^x{x — z). 

Anonfac = ^/fi/D/TrVTas (yU) (Ci + — ^) — J / dx dy dz (f)B{x)(f)D{y) 

V2 JVC I\c Jo Jo Jo 



((x - y)m% + 'mcmD)(p{z) + i^^niBix - z) ( 0^(2) 



+ { {-2x + y + z)mB(t){z) - fi„{x - z){ (t)^{z) 



6 



ml 



6 



m 



B 



Duk^ 



(22) 



where Dc = —m\xx — ml and = m\{x — y){x — z). For the annihilation amplitude Aanni, 
it is the same as the one in Eq. (fT9|) since the two modes D^n^ and D^n^ have the same 
annihilation topological diagrams. 
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For the doubly Cabibbo-suppressed decay mode B — > D 7r+, its decay amplitude can be 
written as 

D-TX+) = VubV^aiAfac + Anonfac + Aanni) , (23) 

here, A^ac^ Anonfac and Aanni cau be obtained from the ones of decay mode D^n^ by 

simply exchanging the Wilson coefficients Ci and C2. 

For the B^ D^tx~ decay, its amplitude can be yielded by using the isospin relation 
A{B~ Z^Oyr-) = A(^ D+tt-) - V^Ai^ D^n^). 



IV. TREATMENTS FOR PHYSICAL-REGION SINGULARITIES 



To perform a numerical calculation of the decay amplitudes of i? — Dir decays, the light- 
cone projectors of mesons are found to be very useful, and the details of these quantities are 
presented in the Appendix. Where one encounters the endpoint divergences stemming from 
the convolution integrals of the meson distribution amplitudes with the hard kernels, which 
is caused by the collinear approximation. To regulate such an infrared divergence, we may 
introduce an intrinsic mass scale realized in the symmetry-preserving loop regularization lll.ll2l]. 
At the tree level, it is equivalent to adopt an effective dynamical gluon mass in the propagator. 
Practically, such a gluon mass scale has been used to regulate the infrared divergences in the 



soft endpoint region [13|, [ij, ll5| 



ln( 



ln(- 



(24) 



The use of this effective gluon propagator is supported by the lattice [16| and the field theoretical 



studies 



17|, which have shown that the gluon propagator is not divergent as fast as p-- Taking 



the hadronic scale A = Aqcd, the dynamical gluon mass scale can be determined from one of 
the well measured decay mode. Numerically, we will see that taking Aqcd = 300 MeV, the 
dynamical gluon mass scale is around fig = (1.5 ± 0. 2) Aqcd- 

Another physical-region singularity arises from the on mass-shell quark propagators. It 
can be easily checked that each Feynman diagram contributing to a given matrix element is 
entirely real unless some denominators vanish with a physical-region singularity, so that the ie 
prescription for treating the poles becomes relevant. In other words, a Feynman diagram will 
yield an imaginary part for the decay amplitudes only when the virtual particles in the diagram 
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become on mass-shell, thus the diagram may be considered as a genuine physical process. The 
Cutkosky rules [l^ give a compact expression for the discontinuity across the cut arising from 
a physical-region singularity. When applying the Cutkosky rules to deal with a physical-region 
singularity of quark propagators, the following formula holds 

1 



P 



{ki - k2 - k^Y + ie 



P 



{ki-k2-k,)\ 
1 



-tnS[{k, - k2 - ksYl (25) 
t7T6[ip,-k2-k,f-ml], (26) 



{Pb - ^2 - k^y -ml + ie [(pt -k2- k^Y - ml_ 

where P denotes the principle- value prescription. The role of the 5 function is to put the parti- 
cles corresponding to the intermediate state on their positive energy mass-shell, so that in the 
physical region, the individual Feynman diagram satisfies the unitarity condition. Equations 
(l25l) and fl26|) will be applied to the quark propagators D^a and D^^a in Equation (fT9|) . respec- 
tively. It is then seen that the possible large imaginary parts arise from the virtual quarks 
across their mass shells as physical-region singularities. 

V. NUMERICAL RESULTS 

It is seen that for theoretical predictions it depends on many input parameters, such as the 

Wilson coefficient functions, the CKM matrix elements, the hadronic parameters, and so on. 

□ 

To carry out a numerical calculation, we take the following input parameters pj| 

Ci = 1.117(1.073), C2 = -0.267(-0.179), tub = 5.28 GeV, mo = 1.87 GeV, 

m^± = 139.6 MeV, m^o = 135 MeV, = 4.66 GeV, = 1.47 GeV, 

fso = 216 ± 19 MeV, fo = 223 ± 17 MeV, U = 130.1 MeV, tbo = 1.536 ps, (27) 



Tb- = 1.638 ps, = 1 - AV2, Kf, = AX^ip - zr/), = -A, 



= AW 

The Wolfenstein parameters of the CKM matrix elements are taken as [20(: A = 0.2272 ± 
0.001, A = 0.806 ± 0.014, p = 0.195lJ5:[J^^, f] = 0.326to°oll with p = p(l - ^), f] = ri{l - ^). 
The coefficient of the twist-3 distribution amplitude of the pseudoscalar vr meson is chosen as 
/i^ = 1.5 ±0.2 GeV Q,!^- 

With the above values for the input parameters, we are able to calculate the contributions 
of different amplitudes for each decay mode. Our final results at mh/2 scale are presented in 
Tabled 
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TABLE I: Numerical results at m;,/2 scale of the amplitudes for different diagrams in i3 — > Dtt decays. 
Amplitudes Afac-, -^nonfac-, and Aanni represent the factorizable ((a)and (b) diagrams in Figs. [T] or [2]), 
the non-factorizable ((c) and (d) diagrams in Figs. [T]or [2|), and the annihilation (diagrams in Fig. [3]) 
contributions, respectively. 



Decay 


modes 


Afac 


A-nonfac 


A 






-2.2655 - O.OOeOi 


0.2613 - 0.1862i 


-0.0047 + 0.0053i 






1.7865 - 0.0556i 


-0.0447 + 0.0970i 


0.0008 - 0.0017i 






-0.0603 + O.OOlli 


0.6991 - 0.0022i 


0.0010 + 0.0009i 






-2.3375 - 0.2310i 


-0.4884 -0.1484i 






As a consequence, we are led to the predictions for the quantities z and 5, as well as the 
branching ratios of all the B Dtx decay modes. We present our "default results" of branch- 
ing ratios and detailed error estimates corresponding to the different theoretical uncertainties 
caused from the above input parameters in Tables [TTl and IIIIl respectively. The errors consist 
of three parts: the first one refers to the variation of the dynamical gluon mass scale; the sec- 
ond one arises from the uncertainty due to the CKM parameters A, A, p, and f/; the third one 
originates from the uncertainties due to the meson decay constants and the parameter /i^r. 
From the numerical results given in Tables [I]- llllt we arrive at the following observations: 
(i) For the color-allowed (also Cabibbo-favored) decay mode 5^ — D+tt", the factorizable 
contribution Afac dominates the total decay amplitude, while the contributions from Anonfac 
and Aanni are small. In particular, the contribution of Aanni is so small that we can safely 
neglect it in this decay mode. With the considered uncertainties, it is seen that our result 
is in agreement with the experimental data 22, l23[, and also consistent with the one given 



m 



2^: B{B^ D+TV~) = {2.74:toij) x 10"^ within the allowed theoretical uncertainties. The 



decay amplitude of its CP-conjugate decay mode D tt^ can be obtained from that of 

B D+Tx^ by changing the CKM element VcbV*^ to V^Vud- Since these CKM elements are 
purely real, the branching ratio of 5° D^tx^ decay is the same as that of B^ D^ti~ decay. 

(ii) For the doubly Cabibbo-suppressed decay mode 5^ D^n^, the contributions from 
Anonfac and Aanni are also much smaller than the one from Afac- As the contributions are all 
proportional to the small CKM elements iKifeKdl ~ branching ratio of this decay mode 
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TABLE II: The branching ratios of i? — > Dtt decays with the default input parameters. The theoretical 
results in the second and the third lines correspond to the predictions at the nib/ 2 and scales, 
respectively. The results correspond to fig = 440 MeV. 



Decay modes 


Br{mb/2) 


Br{mb) 


Experiment 


Ref. [24] 


B° ^ Z)+7r-(10-3) 


2.67 


2.20 


2.68 ± 0.12 ± 0.24 ± 0.12 [22j 
2.63 ± 0.05 ± 0.22 [23] 


9 7/1+0.39 
^•'^-0.37 


B- D°TT-{W-'^) 


4.87 


4.73 


4.97 ± 0.12 ± 0.29 ± 0.22 [22] 
4.90 ±0.07 ±0.23 [23] 


.0+0.48 
"3-^"5-0.47 


^ Z)07r0(10-4) 


2.17 


2.04 


2.89 ±0.29 ±0.38 [25] 
2.25 ±0.14 ±0.35[26] 


2.5 ±0.1 




5.54 


4.78 







TABLE III: The branching ratios of i? ^ Dn decays. The theoretical errors shown from left to right 
correspond to the uncertainties referred to as "dynamical gluon mass scale(upper one corresponding 
to fig = 420 MeV and the below one fig = 460 MeV)", "CKM parameters", and "decay constants and 
the parameter specified in the text. 



Decay modes Br (nib /2) Br (nib) 



^- 


-> D + TT- 


-(10 




R7+0.47 +0.62 +0.91 
^•"J' -0.36-0.57-0.85 


9 9^+0.52+0.16+0.80 
^•^"-0.51-0.42-0.62 


B- - 




-(10 




A 07 +0.52 +0.84 +1.74 
-0.48-0.50-0.31 


A 70 +0.63 +0.34 +1.29 
"i-'-J -0.44 -0.89 -1.18 






(10- 


') 


9 17+0.53+0.46+1.14 
-0.45-0.46-0.89 


9 n/l +0.56+0.19+0.95 
^■^^-0.33-0.50-0.75 




D-TT- 


^(10 




cr +0.68 +1.02 +1.52 
'J-'J'^ -0.44 -0.97 -1.38 


A 70 +0.62 +0.81 +1.76 
^•'°-0.48-1.10-1.38 



is found to be at the order of 10"^, and much smaller than that of the Cabibbo-favored decay 
mode. Since the imaginary part of the dominated amplitude Afac is much smaller than the 
real part, the branching ratio of the CP-conjugate decay mode D^7i~ is approximately 

equal to that of the b'^ D^ir^ decay. 

(iii) For the color-suppressed decay mode 'B^ D^n^, the contribution from Anonfac dom- 
inates the total decay amplitude. The result at the scale is smaller than that at the 
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mb/2 scale, but both are in agreement with the prediction given in 2j]: B{B^ D^n^^ 



(2.5 ±0.1) X 10"^ T 
tal data reported in 



le present central value at the mb/2 scale agree well with the experimen- 



26| . but slightly smaller than the recent experimental data given in 25|. 



While when considering their respective uncertainties, our prediction is still consistent with the 
experimental data. 

(iv) For the B~ D^n^ decay, the main contribution originates from the factorizable one 
Afac- Although its decay amplitude can be written as the color-favored — > D^7i~ minus 
the color suppressed B^ D'^n'^ decays, the branching ratio of this decay mode is enhanced 
compared to that of the latter two. The central values of our prediction are well consistent 



with the experimental data given in 22, l23| . On the other hand, when taking into account of 



the theoretical uncertainties, our prediction is also consistent with the one given by the pQCD 
method y : B{B~ D^tt") = (5.43t[]i?) x lO'^. 

(v) Although the branching ratios at the rrii, scale are smaller than those at the mb/2 scale 
in all these decay modes, we can see that the final results have only a marginal dependence 
on the renormalization scale. As for the theoretical uncertainties in these decay modes, the 
errors originating from the dynamical gluon mass scale fig, the CKM matrix elements are 
comparable with each other when fig G (420,460) MeV. However, the uncertainty originating 
from the decay constants and /i^ are dominate in these decays, especially in B~ D^ti" and 
B^ D~TT^ modes. 

It is also interesting to note that the B ^ D transition form factor, F^^^ = 0.634, extracted 
from (a) and (b) diagrams in Fig. [1] is in good agreement with the ones obtained from the other 
methods, such as: F^^^(O) = 0.648 (pQCD method |27[); F^^^(O) = 0.690 (Bauer-Stech- 
Wirbel(BSW) model F^^^(O) = 0.636 (Neubert-Stech(NS) model 28]). 



We now turn to discuss the CP asymmetries in B ^ Dtt decays. As has already been 
discussed above, there are no direct CP violations in all these decay modes. In the following 
discussions, we focus mainly on the time-dependent CP asymmetries of -B ^ D^n^ decays. 

Using the relevant formulas presented in the previous sections, we can predict the CP 
asymmetries in B ^ D^n^ decays and constrain the CKM angle 2/3 + 7 through the two 
observables a and c. Firstly, we present the results of the quantities z and 6 in Table [TVi Taking 
the current constraints for the weak angles (3 and 7 in the SM, we present our predictions for 
the CP asymmetries So+tt~ and S£)-tt+, as well as the two observables a and c. Secondly, 
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taking the weak angles (3 and 7 as free parameters, we show the dependence of the parameters 
Sd+tt- 1 Sr)-n+ and the observables a, c on the angle 2/3 + 7 in Figs. [Hand [5], respectively. 

TABLE IV: The CP asymmetries for B D^tt^ decays. The results in the middle row denote 
our theoretical predictions for each quantity. The center values correspond to fig = I.SAqcdj and 
the error bars originate from the dynamical gluon mass scale with the upper one corresponding to 
Hg = I.SAqcd and the below one fig = 1.7Aqc£). The results in last row are the experimental data. 



results 


z 5 


Sd+tt- 




a 


c 


Theor 
Exp [29] 








-n ni 7+0-003 
u-'-'J^' -0.001 

-0.030 ± 0.017 


-0.022 ± 0.021 



From Table IIVI and Figs. H] and we come to the following observations. 

(i) For B D'^'K^ decay modes, although there are large strong phase difference between 
the Cabibbo-suppressed and the Cabibbo-favored decay amplitudes, the CP- violating 
parameters 3^+^^- and S'£)-^+ are found to be small (—0.01 ~ —0.02) due to the smallness 
of the ratio z ( ~ 0.02). In addition, due to our predictions for the two parameters Sd+t^- 
and SD-n+ are comparable to each other for a given dynamical gluon mass scale, the 
value for the parameter c is nearly zero. 

(ii) The CP- violating parameters Sd+tt- and Sd-tt+ are not sensitive to the choice of the 
dynamical gluon mass scale, especially when the dynamical gluon mass scale is chosen 
above the central value I.SAqcd- However, both of them have a strong dependence on 
the weak angle 2(3 + 7. The same conclusion is also applied to the two observables a 
and c. Our predictions for the two observables are consistent with experimental data 
when considering the corresponding uncertainties. Unfortunately, it is found that with 
the angle 2j3 + '~f varying within the range (0, 180°), almost all of the values for a and c are 
in the range of the experimental data, which indicates that although direct constraints on 
the angle 2/3-1-7 could be obtained through these parameters, the present experimental 
accuracy is insufficient to improve the knowledge of the apex in the unitarity plane. It is 
expected that more precise measurements in future experiments allow us to extract the 
angle 2(3 + 7. 
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FIG. 5: The same as Fig. HI but for CP observables a and c. 



VI. CONCLUSIONS 



In summary, we have calculated the decay amplitudes, strong phases, branching ratios, 
and CP asymmetries for the B Dtt decays, including both the color-allowed and the color- 
suppressed modes. It has been shown that these decay modes are theoretically clean as there are 
no penguin contributions. As a consequence, direct CP violations are absent. The contributions 
from the factorizable diagrams dominate all the decay amplitudes except for the 5^ D^-k^ 
process. All our predictions for branching ratios are consistent with the existing measurements. 
For the D^n^ mode, our predictions will be faced with the future experiments as no 

data are available at present. Due to small interference effects between the Cabibbo-suppressed 
and the Cabibbo-favored amplitudes, the non-zero CP-violating parameters Sd+t^- and So-n+ 
have been predicted in the B D'^n^ decay modes. It has been shown that the CP- violating 
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parameters have a strong dependence on the weak phase 2/5+7, but they are not sensitive to the 
dynamical gluon mass scale. With the angle 2/3 + 7 varying within the range (0, 180°), almost 
all of the values for the CP-violating parameters a and c are within the range of the current 
experimental data. Thus no constraints on the weak phase 2/3 + 7 could be obtained through 
those parameters based on the current experiment data, and more precise measurements are 
needed in future experiments. 

In this paper, we have further shown that the divergence treatments used in our previous 
work 1^ are reliable. Namely, the endpoint divergence caused by the soft collinear approx- 
imation in gluon propagator could be simply avoided by adopting the Cornwall prescription 
for the gluon propagator with a dynamical mass scale. Note that when the intrinsic mass 
is appropriately introduced, it may not spoil the gauge symmetry as shown recently in the 
symmetry-preserving loop regularization ll| . Meanwhile, for the physical- region singularity of 
the on-mass-shell quark propagators, it can well be treated by applying for the Cutkosky rules. 
The combination of these two treatments for the endpoint divergences of gluon propagator and 
the physical-region singularity of the quark propagators enables us to obtain reasonable re- 
sults, which are consistent with the existing experimental data and also in agreement with the 



ones 



24| obtained by using the pQCD approach. However, this is different from the treatment 
of the latter, where fc^ and Sudakov factors have been used to avoid the endpoint divergence. 

It is noted that the resulting predictions for the branching ratios are in general scale depen- 
dence on the dynamical gluon mass which plays the role of the IR cut-off. This dependence 
should in principle be compensated from the possible scale in the wave functions which charac- 
terizes the nonperturbative effects. In our approach, the dynamical gluon mass may be regarded 
as a universal scale to be fixed from one of the decay modes. For instance, in our present con- 
siderations, if the decay mode D~^it~ is taken to extract the dynamical gluon mass scale, 
we have fig ~ 440 MeV, and the resulting predictions for other decay modes can serve as a 
consistent check. Within the current experimental errors and theoretical uncertainties for some 
relevant parameters, it is seen that our treatment is reliable. In order to further check the 
validity of the gluon-mass regulator method adopted to deal with the endpoint divergence, it is 
useful to extend this method to more decay modes. Anyway, the treatments presented in this 
paper may enhance its predictive power for analyzing non-leptonic B-meson decays. 
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Appendix: Detail calculations of the B Dtt decays amplitudes 

To evaluate the hadronic matrix elements of -B ^ Dtt decays, the meson light-cone distri- 
bution amplitudes play an important role. In the heavy quark limit, the light-cone projectors 
for B, D and vr mesons in momentum space can be expressed, respectively, as 



-ML = ^ { A 75 m - /ip75 (^M^) - ^cr,.n%-^ + ^a,^P,^ '^'''^^ ^^^^ J |^ , (28) 

From the Feynman diagrams shown in Figs. [T]l3l we can get the amplitudes for each decay mode 
using the relevant Feynman rules and the light-cone projectors listed in Eqs. (!28l) . 

For the tree diagrams of b'^ D^n^ mode shown in Fig. [1], the amplitudes of each diagrams 
can be written as 



b ITT'b 



= -^fn9!^j^^r[M''rM'' /'3(l-75)(^5 +m,)7j 
Nc Dbk^ 



= -lUgl^^TrlM^^'-M^iMc +me) A (1-75)] 

fd 

= -^?'^:D^Tr[M-7" ^rf7ni-75)]Tr[A^^7«-M%(l-75 
= -^'^:D^Tr[Al'^7^(l -75) ^„7lTr[A^^«A^^7.(l -75)], (29) 
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where stands for the jth(j = a, b, c, d) diagrams in FigH], km and k the momentum of m 
quark propagator and gluon propagator, respectively. Furthermore, Dm and k^ represent for 
the m quark propagator and gluon propagator, respectively. 

In Fig. [T](a), the vr meson can be written as a decay constant since it originates from the 
vacuum. Inversing the fermi lines and writing down the B meson projector Ai^, gluon vertex 
—igs'J°'T^j, D meson projector M.^ ^ the four quark vertex 7^(1 — 75), b quark propagator 



and another gluon vertex —igs'j^T^i in a trace one by one, and finally the gluon propagator 
-'9c^f-b ^ ^j^g amplitude M". y# can be calculated in a similar way. In Fig. [I](c), the 

TT meson can no longer be written as a decay constant any more since it exchanges a gluon with 



the spectator quark. Writing down the vr meson projector A4^, gluon vertex —igs'~f°'T^-, d quark 
propagator ^ and the four quark vertex 7^(1 — 75) in turn in one trace, and writing down the 
B meson projector M.^ , gluon vertex —igs'y^T^i, D meson projector M.^ and the four quark 
vertex 7^(1 — 75) in the other trace one by one, and finally the gluon propagator we 
can get the amplitude Similarly, we can get the amplitude Summing up the former 
and the latter two quantities in Eq. (129|) . we can get the factorizable part Afac (Eq. (fT4l) ) and 
the nonfactorizable Anonfac (Eq. ffTSl) ). respectively. 

As for the annihilation diagrams for 5° D^iv^ in Fig. [HI the amplitudes can be written 



as 



^ = z/BPrTr[A^^(-z^.7"7-.)Tr^7.(l-75)A<-(-^^?.7^T,^j]^^^'^^ 

.C_F 1 

= ^/^PfTr[^%(l-75)^(-^^?s7"7;^)A^'^(-^^?s7'^T,^J]^^^. 



I\r< JJrn.Kr, 



-lf^gl^-l-i:r[M'> A (1 -75) ^«a7"-M-7a], 

T.[M^{-rg,rT^^^^^{l-^,)]Tr[M^^,{l-^,)^^ 
-^?,^^-^Tr [A^ V ^.a7^(l - 75)] Tr [^^7,(1 - 75)-M'^7j , 



^ = Tr[A<V(l-75)Tr4 (-^^?^7"^1)]Tr[-M%(l-75)-M'^(-z^?.7''T^^l^^ 



-^?,^^-iT^Tr[A^V(l -75)(^.a +m,)7"]Tr[^%(l -75)-M^7j, (30) 



where kma and ka stand for the momentum of m quark propagator and gluon propagator, and 
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Dma and k'^ represent for the m quark propagator and gluon propagator in these annihilation 
diagrams, respectively. Summing up the four quantities in Eq. (!30|) . we can get the annihilation 
contribution Aanni (Eq- f|T9l) ) of this decay mode. 

Similarly, as for the tree diagrams of D^n^ decay mode in Fig [21 its amplitudes can 

be written as follows 



= -^/z)^?,'^-^Tr[A^VA^" A(l -75)(^6 +m,)7„], 
Nc Dbk^ 



= -2/z,^?,'-^y^Tr[A^V-M'^7a h A (1-75 
Nc Ddk^ 

M = Tr[M''Hg,rTt,)jrr^ .7'^(1 - 75)]Tr [A^^(-z^.7'^T,^,)A<-7.(l - 75)] 

= -^^?'^:D^Tr[AlV(^c +m,)7''(l-75)]Tr[AlVA^-7M(l-75 



k^ 



= -^^7,^^-^Tr[A^V(l-75) ^.7lTr[A<^7.-M'^7;.(l-75)]. (31) 

We can get the factorizable contribution Afac (Eq. ( |2T]) ) and the nonfactorizable part 
A-nonfac (Eq. (l22l) ) by summing up the formerand the latter two quantities in Eq. (I3T1) . 

As for the annihilation diagrams for 'B^ D^ir^ decay, its amplitude is the same as the 
one in Eq. (!30!) since the two modes D^-k^ and D'^tx" have the same annihilation topological 
diagrams, which are shown in Fig. [3l 

For the doubly Cabibbo-suppressed decay mode 5^ D^ir'^, its decay amplitude can be 
similarly expressed as the ones in Eq. ( 13T1) due to the same topological structure in these two 
decay modes. 

Finally, for the B~ D^ir^ decay, since its Feynman diagrams are the sums of the color- 
allowed and the color-suppressed one, we can easily get its amplitudes using the above results. 
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